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We present a brief review of current methods for the calculation of multi-loop amplitudes including recent develop- 
ments. As an example we present the calculation of the second moment of the heavy quark current correlator and the 
extraction of the values of the charm and bottom quark masses using the results of this calculation. 



1. INTRODUCTION AND MOTIVATION 



In the era of the LHC and in preparation for the ILC there is need for multi-loop calculations to match the 
precision of these experiments from the theory side. The classes of problems currently accessible reach from multi-leg 
multi-scale problems at one-loop up to no-scale propagator, and tadpole diagrams at four-loop level. Many problems 
fall into these categories directly while others can be cast into the required form by appropriate expansions in small 
parameters. A typical multi-loop calculation proceeds as follows: The needed diagrams are generated and if necessary 
expanded to decrease the difficulty of the problem. The resulting Feynman integrals are then reduced to a small 
set of simpler integrals, which have to be calculated in a further step. While the reduction of the integrals can be 
automatized completely and is mostly independent of the problem at hand, the calculation of the remaining integrals 
must still be done manually for every calculation. In the following we give a brief review of the current status of the 
reduction to master integrals which with recent developments is now in a good shape such that the calculation of the 
master integrals now again poses the most difficulties in multi-loop calculations. 



One of the main ingredients of many multi-loop calculations are the so-called Intcgration-By-Parts (IBP) identities 



which hold in dimensional regularization. Here we consider typical Feynman integrals, where Zj are loop and Pi are 
external momenta. These identities can be used to relate different integrals. Using these relations the full set of 
needed integrals, which can be of the order of 10 4 — 10 6 integrals, can be reduced to a small set of so-called master 
integrals. In a typical calculation there appear O(10 — 100) master integrals, which remain to be calculated directly. 
Another main advantage of the reduction to master integrals is that many properties like gauge invariance can be 
checked on the level of master integrals, since these are linearly independent. 

To make use of the IBP identities the Laporta algorithm 0], a Gauss-elimination-like algorithm, is frequently 
used. In order to use the Laporta algorithm, one has to choose an ordering to ensure that more difficult integrals 
are expressed through simpler ones. Once this ordering has been chosen the algorithm proceeds as follows: Choose 
a set of integrals, which you like to solve; pick an integral from the set and generate an IBP; solve for the most 
difficult integral in the equation and continue with the next IBP or integral. The system of equations will become 
overdetermined if one chooses a sufficiently large set of integrals, since there are Nl(Nl + 2Ne + 1) /2 IBP identities 
generated by a given integral, where Nl and Ne are the number of loops and legs, respectively. 

Although this method has been successfully applied in many calculations it has some drawbacks. The complexity 
of the Laporta algorithm is 0(n 3 ), where n is the number of equations. In combination with the combinatorial 
properties of the problem this leads to fast growth of the size of the problem to be solved. Furthermore in typical 



2. TECHNICALITIES 
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applications the system of equations is overdetermined by a factor 3-5, so quite some time is spent just checking that 
the system of equations is consistent. 

A different approach for solving the set of IBP identities is provided by Groebner bases. Groebner bases arose in 
the context of ideals of rings over multivariate polynomials. More precisely, they were introduced while trying to 
solve the problem, how to decide if a certain polynomial belongs to an ideal of a ring of polynomials. An ideal I is a 
subset of the ring R such that for every q G I,r G R, the product is an element of /, q ■ r G /. 

An ideal can be generated by a set of polynomials Q = {qi, . . . , q n } using I = {p G R\p = ^ r i1ii r i £ R\- Given 
these generating polynomials, the question whether a given polynomial p G R is an element of the ideal, can easily 
be answered in the case of univariate polynomials, where the operation of polynomial division is unambiguous. If 
the polynomial division by the generating polynomials leaves no remainder, p is an element of the ideal. 

This statement docs no longer hold, when dealing with multivariate polynomials. Here the remainder of multivariate 
divisions by multiple polynomials depends on the chosen ordering of the polynomials. Therefore if the remainder of 
such divisions is not zero one may not conclude that the polynomial is not an element of the ideal. 

At this point Groebner bases prove very useful since a multivariate division with respect to elements of a Groebner 
basis is again unambiguous. This leads to one of the definitions of a Groebner basis: A set B C R is a Groebner basis 
if the reduction, i.e. multivariate division, of any r G R with respect to the elements of B yields the same remainder 
independent of the chosen order. Using the Buchberger algorithm a Groebner basis can be constructed from any set 
of polynomials generating the ideal. 

This knowledge can be used for the solution of the IBP-identities. For this purpose the IBP identities have to be 
written in terms of shift operators 0, 0, or differential operators 0, 0| ■ The IBP relations then form a left-ideal in 
the now non-commutative ring of shift operators for which a Groebner basis has to be constructed. Every integral 
can be obtained by applying shift operators to a basis integral, i.e., J{a,b) — pj(l,l),p G R. Since every element 
p of the ring can be expressed in the form p = r iPi + 9j where pi G / , 7"j, q G R the integral can be reduced via 
J(a, b) = pJ(l, 1) = (X) r iPi + <?)>/(l, 1) = qJ(l, 1), since riPi G / and therefore riPiJ(l, 1) = 0. 

Although this procedure does, in principle, generate a reduction of any integral there are still some problems. The 
Buchberger algorithm is guaranteed to generate a basis but it might need a very long time and a lot of resources. 
In practice there is no implementation of the Buchberger algorithm which succeeds in constructing a basis for non- 
trivial problems. Furthermore, as can be checked in simple cases, the application of Groebner bases to the reduction 
problem does not yield the minimal number of master integrals that can be reached using a Laporta-like kind of 
reduction method. 

To circumvent the problematic issues of the Groebner basis approach, a modification of the Buchberger algorithm 
has been introduced to construct so-called s-bases. Unfortunately in this approach the construction of a basis is no 
longer guaranteed to succeed and is highly dependent on the chosen ordering. For certain topologies a basis cannot 
be found and therefore the method cannot be applied. In this cases a conventional Laporta approach still has to be 
used in the sectors where no basis can be found. 



3. EXAMPLE AND APPLICATION 

As an example for the power of the methods outlined in the previous Section we present the calculation of the 
second moment of the vector correlator at 0{a i s ) and its application to the determination of the charm and bottom 
quark masses from experimental data. 

The charm and bottom masses can be extracted from the i?-ratio R(s) = - 1 - ° T^-^ ^ The extraction 

of the quark masses can be done by considering the moments C^ xp of the ratio defined by 

C n P = f^ds. (2) 
These moments are related to the Taylor coefficients C^ 160 of the low energy expansion of the heavy quark current 
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correlator 



(jLX 

\AM 2 J 



where II is the vacuum polarization function of the photon 

n r («V" ~ <fq v ) = IP" 



n(<z 2 ), 



(0\Tf( q )f(0)\0) 



(3) 



(4) 



with the heavy quark current j 11 = tpj^ip. These theoretical moments can be calculated in an expansion in the strong 
coupling constant a s in perturbative QCD 
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theo 
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To perform the calculation the initially present propagator- type diagrams are expanded in q 2 /AM 2 where q is the 
external photon momentum and M the mass of the heavy quark. This expansion leads to tadpole diagrams to which 
the techniques presented in the previous Section can be applied. 

.30- 

At four loops the 



The Taylor coefficients have been calculated at 0{a 2 ), i.e. three-loop order, in 



first moment has been calculated in 



1( 
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and the second moment in [17| . see [1 



19j for higher moments for 
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211 ] . can be 



purely fermionic contributions. Using the second moment the analysis presented in [201 ] . see also 
upgraded to include the now complete knowledge of the second moment. Using this newly available information the 
error of the bottom quark mass can be reduced to obtain 



m b (10GeV) = 3.607(19)GeV, 
which has to be compared with old value of 

TO fc (10GeV) = 3.609(25)GeV. 

In the case of the charm quark one finds 

m c (3GeV) = 0.976(16)GeV 

instead of 

m c (3GeV) = 0.979(22)GeV. 



(6) 



(7) 



(8) 



(9) 



Even though the shift in the absolute value is quite small the error gets reduced by around 25%. 

Since the different moment weights the experimental input differently precise knowledge of the second moment 
allows for a consistency check of the method used for the extraction of the quark masses. Especially in the case of 
the bottom quark the second moment is favored over the first one since it is more sensitive to the threshold region 
and the narrow resonances. 



4. CONCLUSION 

We presented a short review of the current status of an integral part of multi-loop calculations, namely the reduction 
to master integrals. Using recent Groebner basis inspired methods in combination with Laporta-like approaches most 
reduction problems can be solved in an efficient way and the most problematic part is again the calculation of the 
master integrals. As an example we presented the calculation of the second moment of the heavy quark current 
correlator at four-loop order and its application to the extraction of the charm- and bottom-quark masses from 
experimental data. 
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